Theoretical uncertainties on quasielastic charged-current neutrino— nucleus cross 

sections 

M. Valverde,-'^ J. E. Amaro,^ and J. Nieves-'^ 

Departamento de Fisica Atomica, Molecular y Nuclear 
Umversidad de Granada, E- 18071 Granada, Spain 

We estimate the theoretical uncertainties of the model developed in Phys. Rev. C70 055503 for 
inclusive quasielastic charged-current neutrino-nucleus reactions at intermediate energies. Besides 
we quantify the deviations of the predictions of this many body framework from those obtained 
within a simple Fermi gas model. A special attention has been paid to the ratio a{fi)/a{e) of 
interest for experiments on atmospheric neutrinos. We show that uncertainties affecting this ratio 
are likely smaller than 5% 
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I. INTRODUCTION 



Recently an interest in neutrino scattering off nuclei has raised because of its implications in the experiments 
on neutrino oscillations based on large Cerenkov detectors. Atmospheric neutrino oscillations are dominantly due 
^^ , to i'^ — > Vr flavor mixing, which has been already confirmed by the Super-Kamiokande experiments [ij. Future 
^^ ' experiments aim at observing sub- leading effects on top of the dominant one g, which are due to atmospheric 
_i I i/e flavor oscillations. Observability of these effects depends strongly on the assumed experimental and theoretical 



systematic uncertainties. The most simple model regarding these effects supposes that these oscillations are controlled 

r»0 [ by neutrino masses (through their mass squared differences Arri^) and by a 3 x 3 unitary mixing matrix similar to 

T^ > the CKM matrix in the quark sector. Thus we are provided, besides the mass squared differences, with 3 (mixing 

^^ ' angles) + 1 (CP-violating phase) parameters not fully determined yet. To reduce the error size on these parameters 

"nI , requires an improvement in already existing neutrino detectors. These experiments use the nuclei in water as target 

for the incoming neutrinos, so a nuclear reaction model is needed for the data analyses. One of the main sources of 

systematical error in these experiments is precisely the cross sections for quasielastic (QE) neutrino scattering J3|. An 

accurate theoretical framework for the neutrino-nucleus dynamics, and a reliable estimate of the uncertainties affecting 

i-G its predictions will definitely help to have a close control over the systematics affecting the oscillation parameters. 

"T"^ There is a rich recent literature on the Charged Current (CC) neutrino-nucleus reactions in the QE region at 

O i' intermediate energies Q-jlJl- However, a systematic analysis of the theoretical uncertainties affecting the predictions 

D of these works is often missing. We will use here the Many Body Framework (MBF) developed in Ref . [Ij| ^ . Starting 

i-C from a Local Fermi Gas (LEG) picture of the nucleus, which automatically accounts for Pauli blocking, several 

^ , nuclear effects are taken into account in that scheme: i) a correct energy balance, using the experimental Q— values, 

is enforced, ii) Coulomb distortion of the charged leptons is implemented by using the so called "modified effective 

momentum approximation", in) medium polarization (RPA), including A— hole degrees of freedom and explicit pion 

j^ • and rho exchanges in the vector-isovector channel of the effective nucleon-nucleon force, and Short Range Correlation 

" " ' (SRC) effects are computed, and finally iv) the nucleon propagators are dressed in the nuclear medium, which amounts 

to work with a LEG of interacting nucleons and it also accounts for reaction mechanisms where the gauge boson, 

W~^ or W~ , is absorbed by two nucleons. This model is a natural extension of previous studies on electron 1^, 

photon [13 and pion [isL Il9l| dynamics in nuclei. Even though the scarce existing CC data involve very low nuclear 

excitation energies, for which specific details of the nuclear structure might play an important role, the model of 

Ref. [iJI provides one of the best existing combined descriptions of the inclusive muon capture in ^^C and of the ^^C 

{v^,H~)X and ^^C (i^e,e~)A" reactions near threshold. Inclusive muon capture from other nuclei is also successfully 

described by the model. Above 80 or 100 MeV of energy transferred to the nucleus, this MBF leads also to excellent 

results for the (e, e') inclusive reaction in nuclei, not only in the QE peak, but also in the A and the dip (located 

between the QE and the A peaks) regions 16]^. It also successfully describes the absorption of real photons by nuclei 

in this energy regime |17|. 



^ This mod el h as been recently extended to the Neutral Current (NC) sector and to the study of nucleon knock-out reactions induced by 

neutrinos Il5l . 
^ Data in ^■^C, '^''Ca and ^"^Pb of differential cross sections for different electron kinematics and split into longitudinal and transverse 

response functions are successfully described. 



In this work we pay special attention to the source and size of the theoretical uncertainties affecting the predictions 
of Ref. [iJI . Firstly we assign to each of the main inputs of the model a reliable uncertainty. For the experimentally 
determined parameters this obviously has to be determined by the experimental error, while for the model dependent 
parameters, we will assume theoretically founded sizes for their errors. Then we propagate the errors by means 
of a numerical simulation, that is: we consider the input parameters to be represented by uncorrelated Gaussian 
distributions, and by means of a Monte Carlo (MC) simulation, we find for any observable predicted by the model its 
derived probability distribution, which specific features will determine its associated theoretical uncertainty. 

Besides, there exist some systematic errors associated to the validity of the hypothesis in which the scheme of 
Ref. [iJI is based. Those are harder to estimate and will be discussed at the end of this paper. 

II. SOURCES OF THEORETICAL ERRORS 

The main inputs of the model of Ref. li| are: 

• Lepton and hadron masses, electro weak coupling constants, which we will assume to be errorless. 

• Neutrino-nucleon form factors. 

The neutrino-nucleon interaction is assumed to be of the V — A type. The vector form factors are related to the 
electromagnetic ones by means of the isospin symmetry, which we will consider exact. For the electromagnetic 
form factors we use the Galster parameterization [2fl|, 



-iN r<N 



(2) 
(3) 



with T = — g^/4Af^, M the nucleon mass, N standing for n (neutron) or p (proton) and /i^v, the corresponding 
nucleon magnetic moment. The axial part is ruled by the axial form factor 
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The pseudoscalar Gp form factor is related to the axial one by means of the partially conserved axial current 
hypothesis (PCAC). The non-pole contribution to Gp, at q^ = 0, is around 200 times smaller 21] than the pole 
contribution taken into account by means of PCAC. We will not explicitly consider this source of systematics, 
though we will assume an uncertainty for gA larger than that generally used in the literature. 

Thus, from these form factors we are left with four input parameters: Mp), A„, gA and Ma, since we neglect 
any type of uncertainty in the nucleon magnetic moments. The Particle Data Group (PDG) compiles several 
determinations of 5yi(0)/.gv'(0)^ from neutron beta decay studies, ranging most of them in the interval 1.25- 
1.27 _22|, we will adopt a conservative point of view and we will take here 1.26±0.01, though the average error 
quoted by the PDG is more than three times smaller. For the axial cutoff mass, we will assume 1.05 ± 0.14 
from the analysis of the vd -^ IJ'~PP reaction carried out in Refs. |23ll24l |. Finally we will take a 5% and 10% 
error on Mp) and A„, respectively. Parameters are compiled in Table Q] 

The nuclear medium effective baryon-baryon interaction used in the computation of the RPA effects. 

The effective nuclear interaction is included in the model starting from an interaction between particle-hole (ph) 
pairs of the Landau-Migdal type 

V ^cq {fo{p)Tif2 + g'a{p)Si32Tif2} , (5) 



^ We assume gv(0) = 1 according to the conserved vector current hypothesis. 



for the isovector channels'*^. The p{r) dependence of /q is hnearly parameterized as [2 

p{r 



• 



/o(pW) = ^/f"^ 



P(0) 



ft^ (6) 



Thus, we have from the tt channel three parameters: cq, /q and Jq ■ In the channel aaff we use an 
interaction with explicit tt (longitudinal) and p (transverse) exchanges and then we replace ,181] 

3 

co9o{p)^i^2TiT2 -> fiT2 ^ crl'^iVfj'' (7) 

where 

Vr/ = mMq) + (% - q,qj)Vt {q)) (8) 

g = q/\q\ is an unitary vector parallel to the transfered momentum and the strengths of the ph-ph interaction 
in the longitudinal and transverse channel are given by 
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Besides, A(1232) degrees of freedom are also taken into account in this channel. The ph-Ah and Ah-Ah 
effective interactions are obtained from the interaction of Eq. {T)) by replacing a ^r S, t ^ T, where S, T 
are the spin, isospin A''A transition operators [ig and f ^ f* — 2.13 /. The SRC functions g'l and g'^ have 
a smooth g— dependence [Ig, which is not considered here, since we will explore only low and intermediate 
energies and momenta, and thus we take g'i{q) — g'til) = g' flJLIlGlll^ . Hence, we end up with six additional 
parameters: /, /* A^r, Cp, Kp and g' . We would like to point out that all these nine parameters, which are 
used to parameterize the medium effective baryon-baryon interaction, were adjusted long time ago [iSl, zg| , and 
since then have been successfully used in several nuclear calculations at intermediate energies [iq- (19] . Here 
we have assumed uncorrelated Gaussian distributions with relative errors of 10%, for all these parameters (see 
Table QJ, except for the constant cq, for which we have not considered any type of uncertainty, because it always 
appears multiplying the parameters /q and /q that have already 10% error within our analysis. 

The nucleon self-energy fS). 

The nucleon self-energy is used to dress the nucleon propagators in the medium. Its real part modifies the 
free nucleon dispersion relation, while the imaginary part takes into account two nucleon absorption reaction 
channels. The result of it is a quenching of the QE peak with respect to the simple ph excitation calculation and 
a spreading of the strength, or widening of the peak. The integrated strength over energies is not much affected 
though. Most of the effect comes from the consideration of the real part of the nucleon selfenergy, as it was first 
pointed out in Ref. [I^^- The model of Ref. (IJ] uses the results of a semiphenomenological approach developed 
by Fernandez de Cordoba and Oset in J2y|. The nucleon selfenergy, 'S{p^,p; p), calculated in this latter reference 
depends explicitly on the nucleon energy and momentum 26] and leads to nucleon spectral functions in good 
agreement with accurate microscopic approaches like the ones of Refs. ^J-.29J. Because in great part the model 
of Ref. [2g is not entirely microscopical, it is hard to identify its parameters. Here we have assumed a 10% 
relative error (Gaussian) for the nucleon selfenergy in the medium. Ten percent error is a reasonable choice, 
since it safely covers the existing differences between the results of Ref. [2g and the microscopic approaches of 
Refs. |27|-|29|. As discussed above, the change in the nucleon dispersion relation is more important than the 
inclusion of the small nucleon width in the medium, related to the quasielastic channels, which will account for 
W^ absorption by two nucleons. Neglecting completely the imaginary part of the nucleon selfenergy leads to a 
considerable reduction in computation time and given the quality of this approximation, it has been used in all 
MC simulations performed to estimate the theoretical uncertainties of the results of Ref. [iJI • 



* The isoscalar channels do not contribute to CC induced reactions. 

^ See for instance solid (red) and dotted (magenta) lines in Fig. 4 of that reference. 



• Proton and Neutron matter densities. 

Charge densities are taken from 31] and proton densities are deduced from them. On the other hand, neutron 
densities are taken approximately equal (but normalized to the number of neutrons) to the proton ones, though 
small changes are considered, inspired by Hartree-Fock calculations with the density-matrix expansion |23| and 
corroborated by pionic atom data '33l. In this work we will present results for oxygen, carbon and argon, for 
the first two we use a modified harmonic oscillator distribution MHO (parameters can be found in Table I of 
Ref. ^^), while for argon we use a two-parameter Fermi distribution (parameters can be found in Table I of 
Ref. [13). We take a 5% relative error (Gaussian uncorrelated) for all parameters, which is about one order 
of magnitude larger than the quoted errors for charge density parameters, but that safely covers uncertainties 
related to the neutron distributions and to the procedure of taking out the the finite size of the proton and 
neutron particles. 



III. RESULTS AND CONCLUDING REMARKS 

By means of a Monte Carlo simulation, we generate a total of 2000 sets of input parameters^ from an uncorrelated 
multidimensional Gaussian distribution, with central values and standard deviations compiled in Table ^ For each 
of the sample sets, we compute the different observables discussed in this work, and thus we obtain the distributions 
of all of them. Theoretical errors and uncertainty bands on the derived quantities are always obtained by discarding 
the highest and lowest 16% of the sample values, to leave a 68% Confidence Level (CL) interval. 

In Figs. Handle we present electron and muon neutrino inclusive QE integrated cross sections from carbon, oxygen 
and argon. Several comments arc in order, 

• As mentioned above, the imaginary part of the nucleon selfenergy has been neglected when performing the 
MC simulation. Let us look at Fig. Q as can be appreciated there, the differences between the central line^ 
of the 68% CL band and the full model prediction, which includes the effects of the imaginary part of E, are 
almost negligible and significantly smaller than the size of the 68% CL errors. This corroborates the findings of 
Ref. | l5| , and though we will make graphically such a comparison in most of the plots that will be presented in 
what follows, we will not make any further comment about it. 

• Nuclear effects beyond implementing Pauli's exclusion principle and enforcing the correct energy balance are 
sizeable and much larger than the uncertainties on the predictions deduced from the MBF presented in Ref. [lj| . 

• In oxygen, we estimate separately the uncertainties (error bars® in the middle plots of Fig.^ due to the imprecise 
knowledge of the nucleon densities and of the parameters entering in the model used in Ref. IHI to compute 
nuclear effects (RPA and nucleon self-energy). The value of the relative uncertainty due to nuclear effects 
decrease with energy, while relative errors induced by the neutrino-nucleon form factors increase with energy, 
and at the higher end, they could be comparable to those affecting the evaluation of the nuclear effects. 

• Uncertainties on the integrated cross sections are of the order of 10-15%, which turn out to be similar to those 
assumed for the input parameters (Table QJ. Hence, predictions of Ref. [Tjl seems stable and no fine tuning 
parameters are identified in the model. 

• As can be appreciated in Fig.|21for two particular cases, theoretical uncertainties on the cross section can be, 
in a good approximation, modeled by a Gaussian distribution. 

Conclusions are similar for differential cross sections. As an example, in Fig. |31 we show the electron neutrino 
inclusive QE double differential cross section, at a fixed lepton momentum transfer, in oxygen as a function of the 
lepton energy transfer. 

Theoretical errors cancel partially out in the ratio cr(/i)/cr(e) = a{iyfj,+^ Z -^ /x~ +X)/a{ve+^ Z -^ e~ +X), as can 
be appreciated in Fig. ^ Theoretical uncertainties on this ratio turn out to be smaller than 1%. On the other hand, 
predictions for this ratio obtained from a simple Lindhard function^ incorporating a correct energy balance in the 



° We have checked that the errors quoted in the following are already stable when 1500 event simulations are performed. 

"^ It has been obtained with the central values of the parameters quoted in Table |l] and neglecting the imaginary part of the nucleon 

selfenergy. 
* Hence, these error bars do not take into account the uncertainties of the neutrino-nucleon form factors. 
^ It is to say from a non-interacting local Fermi gas model of non-relativistic nucleons. 
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FIG. 1: (Color online) Electron and muon neutrino inclusive QE integrated cross sections from carbon, oxygen and argon, as a function 
of the neutrino energy. In all cases non- rela tivistic nucleon kinematics has been employed. Results denoted as "Full model" are obtained 
from the full model developed in Ref. |14| . while those denoted as "Pauli" have been obtained without including RPA, Coulomb and 
nucleon selfenergy effects. We also give the 68% CL band (red or solid lines). For oxygen, the error bars (denoted as "Nuclear") account 
for the uncertainties due to the imprecise knowledge of the nucleon densities and of the parameters entering in the model used (Ref. Il4l ) 
to compute nuclear effects (RPA and nucleon self-energy). 
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FIG. 2: (Color online) MC uncertainty distributions for muon neutrino inclusive QE integrated cross section from argon, at 250 (left) and 
450 (right) MeV incoming neutrino energies. The dashed Hnes stand for Gaussian distributions with central values and variances indicated 
in each panel. 
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FIG. 3: (Color online) Electron neutrino inclusive QE differential cross section, at a fixed lepton momentum transfer of 450 MeV, in 
oxygen as a function of the lepton energy transfer. The incoming neutrino energy is 400 MeV and non-relativistic nucleon kinematics has 
been employed. Results denoted as "Full model" are obtained from the full model developed in Ref. J^, while those denoted as "Pauli" 
have been obtained without including RPA, Coulomb and nucleon selfenergy effects. We also give the 68% CL band (red or solid lines). 



reaction (lines denoted as "Pauli" in the plots) differ from the results obtained from the model of Ref. [lj| only at the 
level of 5%, in sharp contrast with the situation found for each of the the individual cross sections a{v^+^Z -^ fi~ +^) 
and a{i/e +^ Z — > e~ + X) (see Fig. ^. This is consistent with the findings of a somewhat similar analysis carried 
out in Refs. |34L ISSL l3a | . Finally, here we also find that the theoretical uncertainties on the ratio cr(/i)/cr(e) can be, in 
a good approximation, modeled by a Gaussian distribution, as can be seen in Fig. [S] 

For antineutrino induced reactions we find a totally parallel scenario to that discussed so far for neutrino-induced 
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FIG. 4: (Color online) Ratio of inclusive QE cross sections a(fi)/a{e) for carbon, oxygen and argon, as a function of the incoming neutrino 
energy. In all cases non-relativistic nucleon kinematics has been employed. Besides, the 68% CL band (red or solid lines), we also give 
results (crosses) from the full model developed in Ref. Il4ll . and from the latter model without including RPA, Coulomb and nucleon 
selfenergy effects (line denoted as "Pauli"). For oxygen, the error bars have the same meaning as in Figs.^and |21 
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Md 


= 0.843 ± 0.042 GeV 


/c 
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Xn 
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Ma 


= 1.05 ± 0.14 GeV 




/ = 1.00 ± 0.10 


9A 


= 1.26 ± 0.01 




/* = 2.13 ± 0.21 

A^ = 1200 ± 120 MeV 

Cp = 2.0 ± 0.2 

Ap = 2500 ± 250 MeV 

g' = 0.63 ± 0.06 



TABLE I: Central values and errors of the model input parameters, besides we have also included 10% uncertainties (relative) 
in both the real part of the nucleon selfenergy and densities (see text for details). 

ones. 

To finish tliis work, we would like to discuss some other systematic errors associated to the validity of the hypothesis 
in which the scheme of Ref. 11^ is based, and not considered yet. In first place, up to this point a non-relativistic 
nucleon kinematics has been employed. This is because the RPA and nucleon selfenergy models used in Ref. [Ij| 
are based in a non-relativistic effective baryon-baryon interaction in the nuclear medium. The use of relativistic 
kinematics for the nucleons leads to moderate reductions of both neutrino and antineutrino cross sections, ranging 
these reductions in the interval 4-9%, at the intermediate energies considered in this work. Such corrections do not 
depend significantly on the considered nucleus |Ij]. In the ratio a{fj,)/a{e), relativistic nucleon kinematics effects are 
quite small, being always smaller than 1% in the whole neutrino energy interval studied in this work, as can be seen 
in Fig. El 

Second, one might think that a LFG description of the nucleus is poor, and that a proper finite nuclei treatment 
is necessar y. For inclusive processes and nuclear excitation energies of around 100 MeV or higher, the findings of 
Refs. [13, [13 and [13 clearly contradict this conclusion. The reason is that in these circumstances one should sum 
up over several nuclear configurations, both in the discrete and in the continuum, and this inclusive sum is almost 
not sensitive to the details of the nuclear wave function^*^ , in sharp contrast to what happens in the case of exclusive 
processes where the final nucleus is left in a determined nuclear level. On the other hand, the LFG description of the 
nucleus allows for an accurate treatment of the dynamics of the elementary processes (interaction of gauge bosons with 
nucleons, nucleon resonances, and mesons, interaction between nucleons or between mesons and nucleons, etc.) which 
occur inside the nuclear medium. Within a finite nuclei scenario, such a treatment becomes hard to implement, and 
often the dynamics is simplified in order to deal with more elaborated nuclear wave functions. This simplification of 
the dynamics cannot lead to a good description of nuclear inclusive electroweak processes at the intermediate energies 
of interest for future neutrino experiments. 

For all of this, it is sound to assume relative errors of about 10-15% on the QE neutrino-nucleus (differential and 
integrated) cross sections predicted by the model of Ref. JJ|, at intermediate energies. Uncertainties on the ratio 
(7 (l-j.) / (7 (e) would be certainly smaller, likely not larger than about 5%, and mostly coming from deficiencies of the 
LFG picture of the nucleus assumed in Ref. [lj| . 



^^ The results of Ref. |l4l| for the inclusive muon capture in nuclei through the whole periodic table, where the capture widths vary from 
about 4x10'' s-l in ^^c to 1300 xlO'' 8"^ in ^^^Pb, and of the LSND measurements of the ^^c {u^,fi-)X and ^^c {ue,e-)X reactions 
near threshold indicate that the predictions of our scheme, for totally integrated inclusive observables, cou ld eve n b e extended to much 
smaller, of the order of 10 or 20 MeV, nuclear excitation energies. In this respect, the works of Refs. l87|l and |38j for inclusive muon 
capture and radiative pion capture in nuclei, respectively, turn out to be quite enlightening. In those works, continuum shell model 
results are compared to those obtained from a LFG model for several nuclei from '^^C to ^"^Pb. The differential decay width shapes 
predicted for the two set of models are substantially different. Shell model distributions present discrete contributions and in the 
continuum appear sharp scattering resonances. Despite the fact that those distinctive features do not appear in the LFG differential 
decay widths, the totally integrated widths (inclusive observable) obtained from both descriptions of the process do not differ in more 
than 5 or 10%. The typical nuclear excitation energies in muon and radiative pion capture in nuclei are small, of the order of 20 MeV, 
and thus one expects that at higher excitation energies, where one should sum up over a larger number of nuclear final states, the LFG 
predictions for inclusive observables would become even more reliable. 
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FIG. 5: (Color online) MC uncertainty distributions for the ratio of inclusive QE cross sections cr{fi)/a(e) from oxygen, at 300 (left) and 
450 (right) MeV incoming neutrino energies. The dashed lines stand for Gaussian distributions with central values and variances indicated 
in each panel. 
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FIG. 6: (Color online) Relative corrections to the ratio a{fi)/cr{e), in argon, due to the use of relativistic nucleon kinematics. We denote 
by krel ('^norel) the ratio a{ii)/cr(e) deduced from a non-interacting local Fermi gas model of relativistic (non-relativistic) nucleons 
improved by the use of the experimental Q— values (see Ref. 143 for some more details). Thus, knorel would be given by the line 
denoted by "Pauli" in the argon panels of Fig. 1^ 
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